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1. INTRODUCTION 

In modern power systems, it is a mandate of utilizing real-time measurements for better control and 
management of the system. The real-time measurements [1] help us to even protect the devices connected to 
the system during fault conditions. With faster measurement readings taken from the system, a rapid response 
can be created avoiding many fatal conditions. This can be achieved using phasor measuring units (PMU) 
devices [2] connected at different locations of the power system. The placement of these devices can also be 
selected as per the requirement of sensitive devices. PMUs are used for multiple applications including 
assessment of power system stability and security and improving grid protection. These PMUs 
telecommunication devices have an issue of interfacing with external signals introducing disturbance and 
creating errors in the measured signals. These measured signals also have disruptions when the grid system is 
induced by any fault on the transmission lines [3]. 

For the analysis of different data measured by the PMU measuring devices an IEEE 14 bus system 
is considered with fault introduced at different locations of the bus system [4]. The signals from the PMU 
devices are compared with nominal operating conditions with bad data created when a fault is introduced in 
the bus system. Figure 1 shows IEEE 14 bus system, a single-line diagram with multiple sources and loads 
connected at different locations. 
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As observed the bus system [4] comprises different modules connected to different buses. Most of 
the buses have loads that are fed through three generators. Two synchronous condensers are also integrated 
into the system for injection of reactive power improving the voltage profile of the buses. These synchronous 
condensers can be replaced with static condensers which are capacitive VAR devices compensating reactive 
power [5]. The data for the modeling of 14 bus system is as shown in Table 1 to Table 4 as transmission line 
data, static condenser data. Load data and generator rating respectively below. 


Table 1. Transmission line data 


Line From To Line Line 
no. bus bus resistance reactance 
1 1 2 0.01958 0.05927 
2 1 5 0.05513 0.23204 
3 2 3 0.05899 0.18797 
4 2 4 0.06322 0.18632 
5 2 5 0.05595 0.18388 
6 3 4 0.05301 0.18103 
7 4 5 0.02435 0.03511 
8 4 7 0 0.15071 
9 4 9 0 0.66718 
10 5 6 0 0.27802 
11 6 11 0.08498 0.1789 
12 6 12 0.14491 0.28881 
13 6 13 0.08015 0.24227 
14 7 8 0 0.17515 
15 7 9 0 0.12521 
16 9 10 0.02191 0.0825 
17 9 14 0.15761 0.32148 
18 10 11 0.08325 0.18297 
19 12 13 0.22072 0.18828 
Figure 1. Single line diagram of IEEE 14 bus system 20 13 14 0.17063 0.35262 
Table 2. Static condenser data 
Condenser bus no. Qmax (MVAR) 
3 40 MVAR 
9 19 MVAR 
Table 3. Load data 
Busno. P(MW) Q(MVAR) Bus no. P(MW) Q(MVAR) 
1 0 0 8 0 0 
2 21.7 12.7 9 29.5 16.6 
3 94.2 19 10 9 5.8 
4 47.8 -3.9 11 3.5 1.8 
5 7.6 1.6 12 6.1 1.6 
6 112 T5 13 13.5 5.8 
7 0 0 14 14.9 5 


Table 4. Generators ratings 


Generator no. P min P max 
1 10.5MW 165 MW 
2 19.8 MW 82 MW 
3 21 MW 45 MW 


The modeling of the bus system is carried out using the above data with given parameters of 
generators, transmission lines, static condenser, and load data. The total system has an apparent power base 
value (Sbase) of 100 MVA and the voltage of the distribution system is considered to be 11 kV. Therefore, 
the base voltage (Vbase) value is 11 kV which makes the base impedance (Zbase) value given as: 


Z = KV ys 
base MV Abase 
112 
= — = 1.21 ohms 0) 
100 
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with the above impedance value, the transmission lines resistance and inductances are calculated. The 
inductance of the line for the given reactance (X) is given as: 


Xline 
Liine = (2) 


here, w is the angular frequency of the grid system taken as 2.7.50=314.15 rad/sec. 

The above IEEE 14 bus system is integrated with PMU [6] at all the given buses with three-phase 
voltages and currents measurements taken. The system is later included with a fault at one of the transmission 
lines and later updated with unified power quality conditioner (UPQC) module [7] connected at weakest 
bus 14. In this paper, section 1 is included with an introduction to the considered test system and modules 
followed by PMU devices modeling in section 2. In further section 3, the UPQC module configuration is 
discussed. In section 4 the simulation results and analysis are carried out with a comparison of nominal 
operating condition data, bad data during a fault, and updated data with the UPQC module is compared. The 
final section 5 is included with the conclusion to this paper with results discussion followed by references 
used in this paper. 


2. PMU MODULE 

A phasor measuring unit (PMU) is used to monitor transmission line data in real-time mode with 
time synchronization for widely located different substations. It gives magnitude as well as phase angle 
which makes it more advantageous over SCADA [8]. All the measurements in an electrical system vary as 
per the changes taking place in the grid modules. All these measurements are set with benchmark-defined 
values during normal and stable operating conditions. Any data that deviates from the path of these 
benchmark set values [9] during different operating conditions is considered bad data. As mentioned 
previously bad data can be created when the system is introduced with any fault or any external factor 
interference. These faults can be either on the transmission lines or at generators or even on the load side of 
the grid system [10]. The dynamic properties of each bus are noted and can be used for the analysis of the 
system during these faults which can be used for protection and security in the future. 

The main feature of PMU is principal component analysis (PCA) which determines different 
components of the signal. These components can be compared to get bad data generated in the signal during 
disturbances. PCA helps to get lower dimension data which comprises maximum data variables increasing 
the resolution of the data needed [11], [12]. Different measurements can be taken from the PMU module 
comprising voltage magnitudes, current magnitudes, frequencies, or powers. In further detailing of these 
parameters multiple measurements like positive sequence component, negative sequence component, and 
zero sequence components can also be achieved [13], [14]. 

Different measurements have their calculation components which generate data for comparison. In 
MATLAB Simulink software multiple blocks are used to record the data achieved while operating the IEEE 14 
bus system. The main block for measuring the magnitudes of bus voltages and currents is the ‘sequence analyzer’ 
taken from the ‘power system’ toolbox. The parameters of the “sequence analyzer’ are shown in Figure 2. 

As seen the fundamental frequency of the IEEE 14 bus system is considered as 50 Hz and we are 
analyzing the fundamental harmonic component which is n=1. All three components positive, negative and 
zero can be analyzed at specific sampling time Ts. The lower the sampling time Ts the higher the data points 
[15] achieved which helps to store more data for analysis [16]. In the ‘sequence analyzer’ block even the 
angle of the voltages and currents can be determined. Along with this block for measuring the fundamental 
frequency of the buses a phase locked loop (PLL) block is taken from the ‘power systems’ toolbox. The 
parameters of the PLL block are shown in Figure 3. 

As observed the minimum frequency limit is set at 45 Hz for the grid system with a fundamental 
frequency of 50 Hz. The PID regulator gains are set as per the required response of the block. For faster 
response, the rate of change of frequency needs to be increased. The cutoff frequency value mitigates the 
disturbances in the measured data. The recorded data of voltage magnitudes, current magnitudes, 
fundamental frequency, and voltage phase angles of an IEEE 14 bus system when operating in normal 
conditions are shown in the graphs. 

— During normal operating condition: Figure 4 to Figure 7 shows voltage magnitude, current magnitude, 
frequency, and phase angle of all 14 buses respectively during the normal operating condition. All these 
graphs are taken with a simulation time of Isec and the stabilized values are displayed on the right of 
the graph. The IEEE 14 bus system is induced with a fault at one of the transmission lines and the data 
is recorded. A fault is introduced between bus 4 and 5 with a 3-phase to ground symmetrical fault for 
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two cycles from 0.5 to 0.52sec of simulation time. Below are the graphical bad data achieved during 


this condition. 


During fault condition: Figure 8 to Figure 11 shows voltage magnitude, current magnitude, frequency, and 
phase angle of all 14 buses respectively during a fault condition. As seen in the time domain graphs, bad data 
is created during the time 0.5 to 0.55sec due to a fault on the transmission line. For improvement in the 
operating conditions during a fault, an external device is connected to the weakest bus. The system is further 
updated with UPQC and the data is compared with normal, bad, and improvised conditions. 


Sequence Analyzer (mask) (link) 


Compute the positive-, negative-, and zero-sequence components 
of a three-phase signal. 


Parameters 


Fundamental frequency (Hz): 


50 


Harmonic n (1=fundamental): 


1 


Initial input [ Mag, Phase (degrees) ]: 


Sequence: Positive Negative Zero s 


[1, 0] 


Sample time: 


Ts 


Figure 2. Sequence analyzer parameters 


Voltages of 14buses in pu 


0 01 02 03 04 05 06 07 08 09 


Time (seconds) 


Figure 4. Voltage magnitudes of all 14 buses 


03 04 05 06 07 08 09 
Time (seconds) 


Figure 6. Frequency of all 14 buses 


Current magnitudes of all buses 


Phase anlge of phase A all 14 buses 


Parameters 
Minimum frequency (Hz): 
45 


Initial inputs [ Phase (degrees), Frequency (Hz) ]: 
[0, 50] i i o E 
Regulator gains [ Kp, Ki, Kd ]: 

[180, 3200, 1] 


Time constant for derivative action (s): 
ies 


Maximum rate of change of frequency (Hz/s): 
12 


Filter cut-off frequency for frequency measurement (Hz): 
25 


Sample time: 
Ts 


~] Enable automatic gain control 


Figure 3. PLL block parameters 
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Figure 7. Phase angle of all 14 buses 


Bad data analysis and detection using PMU with UPQC integration to ... (Preeti Kabra) 


260 o ISSN: 2088-8694 


n 
2 
2 0.8 3 
= T 
S oe 5 
> N 
= 3 
= z 
© = 
3 °4 = 
s € 
S 2 S 
È 
> 
rs) 
o 
L 4 L L i i 4 i L 0 0.2 04 06 08 1 
o o1 02 03 04 05 06 07 08 09 1 Time (seconds) 
Time (seconds) 
Figure 8. Voltage magnitudes of all 14 buses Figure 9. Current magnitudes of all 14 buses 
a et 
pel 
2 
50.05 2 
S = 
3 o” pE 
v <. 
= 49.95 3 
© g -20 
> a 
5 49.9 5-25 
5 ke] 
fog > 
D 49.85 g °° 
9 -35 
49.8 E 
Q. -40 
49.75 t =t i i i i = + i—i 4 i L k 4 
0 0.1 0.2 03 04 05 06 07 08 09 1 0 02 04 06 08 1 
Time (seconds) Time (seconds) 
Figure 10. Frequency of all 14 buses Figure 11. Phase angle of all 14 buses 


3. UPQC CONFIGURATION 

UPQC [17] is a traditional device with two modules connected to the grid in series and parallel 
connections. The series module is considered a dynamic voltage restorer (DVR) [18] and the shunt module is 
considered a synchronous static compensator (STATCOM). Each module comprises one voltage source 
converter (VSC) which converts AC-DC and DC-AC as per the requirement [19]. Both the VSCs are 
connected to the grid through inductive filters for harmonics reduction. Each VSC comprises 6-IGBT 
switches modeled in 3-legged format for connection to 3-ph lines. The VSC has two sides, a 3-ph AC side, 
and a DC side. The 3-ph AC side is connected to the grid and the DC side is connected to a common DC link 
capacitance for energy storage. The DC link capacitance is a support element for the VSC which helps to 
mitigate different issues in the grid. The DVR module (series) compensates for voltage fluctuations and the 
STATCOM module (shunt) compensates for harmonics [20]. The internal modeling of the UPQC device can 
be seen in Figure 12. 

The 6-switches in each module need to be controlled for the compensation of voltages and 
harmonics in the grid. Both the VSC modules are controlled with synchronized pulse generation concerning 
the source voltages. The DVR controller [21] generates references as voltage components and the 
STATCOM controller generates references as current components. Figure 13 is the controller for the DVR 
module with the Sin PWM technique [22]. 

As observed in Figure 13 the controller takes the input of source voltages (VSabc) and load voltages 
(VLabc) for the generation of series reference voltage components (VSEabc*). To generate these components 
both VSabe and VLabc signals are converted to dq-components [23] and compared to generate the error 
voltages. The Park’s transformation is given as: 
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Va Sin@ -—Cos@ 0][V% 
V% |= |Cos0 Sino OflY%, (3) 
Vo 0 0 LILY 
the series reference dq-components (VSEdq) are given as: 
Vsea = (Kp + SK) (Via — Vsa) — Via — Vsa) (4) 
Vseq = (Kp + J Ki) (Vig =Vsq) = Vig = Vsa) (5) 


here, Kp and Ki are the gain of the PI controller which are set by tuning the module until stabilized voltages 
are achieved. The series reference dq-components are converted to abc components using inverse Park’s 
transformation for Sin PWM technique pulse generation. 


y Sin 0 Cos 0 1 7 
a . 2T 2T d 
p = |sin(0- F) cos(0-7) l k ©) 
V; ; 20 2n V 
c Sin(0 +=) Cos (0 +) 1} Lo 
Compensation Voltage ( V srabe ) 
e E, > Compensation Current(I,, abe) 
oe Load side 
side 


s Q 
ai 5 
Vse k s XS PWM Voltage 
Controller 


Vig’ z Gate Pulses 
a Vig =0 Series VSC 


Figure 13. DVR module control structure 


Similar to this the STATCOM module also has many components taken as input for generation 
reference current components. The controller of the STATCOM [24] is updated with the virtual inertia 
concept which makes this static device operate as a dynamic synchronous condenser for better results 
generation. The updated virtual inertia controller of the STATCOM can be seen in Figure 14. 
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Virtual 


Impedance 


Figure 14. Virtual impedance integrated STATCOM controller 


In Figure 14, controller Vdc* is the reference [25] DC link voltage and Vdc is the measured DC 
link voltage. The Pac reference is given as: 


Pac = (Kp + J K) (Vac — Vac) 7) 
The change in angular frequency is generated concerning Pac and Pac values. 
1 
Aw = (—) (Pac — Pac) (8) 


The reference phase angle for the generation of reference voltage signals is given as: 
8 = f (wn + Awn) (9) 


The magnitude for the reference voltage signals is taken as: 


E* = (K, + J Ki) * (Vise 


-= Vecel-) (10) 


Therefore, the reference current i* components is given as (11). 


r= (e-a ay 


In the above equation e=E*Sin@. L and R are the inductance and resistance of the series inductor connected 
to the shunt VSC [23]. The reference current (i*) is converted to dq-components which are compared to 
measured current igg-components of the VSC output. 

The final reference voltage dq-component is given as (12). 


Vaq = (Kp + J Kr) (ag — tag) (12) 


This final reference voltage dq-component is converted to abc using (6) and the final reference voltage 
signals are fed to Sin PWM technique for controlling the shunt VSC. The IEEE 14 bus system is updated 
with UPQC connected at the weakest bus (bus no. 14) and the results are compared with normal data, bad 
data during fault, and improvised data with UPQC. 


4. SIMULATION RESULTS 

The modeling of the IEEE 14 bus system with UPQC module connected at bus 14 and fault between 
bus 4 and 5 is done with MATLAB Simulink software. The model is run for | sec with fault set from 0.5 to 
0.52 sec in all three models. Figure 15 shows voltage and current magnitudes of all 14 buses during fault 
condition with UPQC and Figure 16 shows frequencies and Phase angles of phase A of all the 14 buses 
during fault conditions with UPQC. Figure 17 shows the comparison of the voltage magnitudes and 
frequency respectively when the data is recorded in normal (blue), fault (red), and integrated with UPQC 
(green) operating conditions for bus 14. It indicates bad data which appeared in normal conditions can be 
improvised with virtual synchronous machine (VSM)-UPQC and ultimately improve the quality of data. 
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Figure 15. Voltages and currents magnitudes of all 14 buses during fault condition with UPQC 
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Figure 16. Frequencies and phase angles of phase A of all the 14 buses during fault conditions with UPQC 
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5. CONCLUSION 


Successful implementation of the IEEE 14 bus system as per the data and integration of the UPQC 
module into the grid system is done. The bus system is introduced with a fault on one of the transmission line 
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creating bad data in the measurements. The bad data is compared to normal data which is achieved during the 
normal operating condition. Four parameters are considered for data recording which include voltage 
magnitudes, current magnitudes, frequencies, and phase angles. All the data generated by the PMU are 
plotted in graphical format for comparison to determine the difference between good and bad data. With 
further updating the bus system with UPQC the new improvised data is also compared showing the 
improvement in voltage magnitude and frequency of the buses. For our analysis, the UPQC device is 
connected to the weakest bus (bus no. 14) and the data achieved is compared. Concluding that the bad data 
created due to fault in the system is stabilized by the UPQC device with lower voltage magnitude drop and 
lower frequency disturbance. The voltage magnitude is maintained near unity and the frequency is remaining 
at 50 Hz even during the fault condition. 
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